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ABSTRACT 



Very low frequency (VLF) electromagnetic wave injection experiments were 
conducted on 23-24 January 1988 from a 42-km horizontal dipole antenna located at Siple 
Station, Antarctica. The experiment consisted of a diagnostic format transmitted for one 
minute every five minutes for a 10 hour period between 1700 UT and 0300 UT. These 
signals were received and recorded at the conjugate magnetic field point location at Lake 
Mistissini, Canada. A detailed analysis of this data clearly demonstrates hot plasma effects 
such as saturated power levels, exponential growth rates, sideband formation and 
triggered emissions due to wave-particle interactions. These hot plasma effects remain 
constant over a time scale of 30 seconds but show large variations over a time scale of 5 
minutes. These VLF signals were used to simulate "whistler waves" which occur naturally 
and are amplified by energetic electrons spiraling around magnetic field lines near the 
geomagnetic equator. Navy VLF communications are strongly affected by the presence of 
whistler waves. 

The electron and whistler wave interaction can be described by a theoretical model 
which is very similar to that used for free electron lasers (FEL). Using computer 
simulation most of the hot plasma effects seen in the Siple Station data can be modeled 
and compared to free electron laser characteristics such as saturation, electron trapping, 
tapering, and sensitivity to energy distributions. 

Two dimensional computer simulations in coordinates z and t have predicted that the 
CEBAF Infrared (IR) FEL can observe limit-cycle behavior when operating within its 
design parameters. The IR FEL is driven by a high quality electron beam with a 
micropulse length comparable to the slippage distance. At moderate values of the 
desynchronism, the optical power will oscillate periodically over several hundred passes 
through the resonator. The limit-cycle power oscillations are caused by "marching 
subpulses" that grow at the trailing edge of the optical pulse through a super-radiant 
process, and pass through the main optical envelope. 
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I. INTRODUCTION 



The study of "whistlers" began with the earliest days of radio communication. 
Whistlers are tones that sweep over many frequencies within the audio range and sound 
like whistles or musical type signals. At the time, prior to World War I, whistlers were 
heard, but their origin was not known. Then in 1919, Barkhausen officially published the 
first account of whistlers while studying the currents in the earth's ionosphere. Between 
1930 and 1935, Barkhausen and Bckersley wrote a paper describing the theory of 
whistling atmospherics in great detail [Ref. 1]. They said, "Whistlers originate from 
impulsive atmospherics, or lightning, and travel through the outer ionosphere, following 
the lines of force of the earth's magnetic field and crossing over the equator at great 
heights." [Ref 1] Even then, they knew that the dispersion measurements of whistlers 
yielded information about the density of electrons in the ionosphere. 

These originally published experimental observations lacked supporting theoretical 
computations, until L. R. O. Storey published his work on whistlers in 1952. Storey 
performed a systematic study of whistlers over many years before making conclusions 
about the variations of whistler properties. Then, he formulated and published the first 
theoretical work on the dispersion relationship, group velocity, and time-delays due to 
field-aligned ducted propagation of a whistler. Finally, in 1965, R. A. Helliwell published 
the first comprehensive book on whistlers that incorporated both experimental 
observations and theoretical explanations for most ionospheric phenomena [Ref. 2]. 

The study of whistlers has long been a topic of concern for the military. Some of the 
first reports of whistlers came from Germans attempting to listen to Allied telephone 
conversions during World War I [Ref. 2]. Since then, whistlers have been a problem to 
communication lines and navigation systems as described in Chapter II. 

In order to better understand whistlers, Siple Station was built in Antarctica in 1972. 
This station was used to transmit very low frequency (VLF) signals along field-aligned 
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magneto spheric ducts which were received at Lake Missiniti, Canada. These VLF signals 
were used to simulate the properties of whistlers and their propagation paths in a 
controlled fashion. Chapter III describes the experiment and equipment used at Siple 
Station. Also included in Chapter III is a brief description of the magnetosphere and the 
Van Allen Belt electron properties. 

Chapter IV begins with a detailed study of frequency and amplitude information 
received from Siple Station on 23-24 January 1988. Examples of single and multi-path 
propagation are shown, and the theory behind ducted propagation is presented. Examples 
of triggered emissions from the transmitted signals are shown and defined. Characteristic 
properties of the wave-particle interaction between the VLF wave and the hot plasma 
electrons are presented. These include exponential growth of the VLF signals, saturation 
of the VLF signals, sideband formation after saturation, and triggered emissions at signal 
termination. The concept of threshold power and the use of Siple Station signals as a hot 
plasma diagnostic tool are presented. 

In Chapter V, a theory is developed to explain the amplification of the VLF signals in 
the inhomogeneous magnetosphere. The whistler mode dispersion relationship is derived 
from Maxwell's equations and the cold electron current. Then the VLF wave equation is 
developed by including the hot electron current in Maxwell's equations. Finally, the self- 
consistent motion of each hot electron is described by solving the Lorentz force equation 
driven by the VLF wave and the earth's inhomogeneous magnetic field. Amplification of 
the VLF signal occurs because of a gyro-resonant condition between the wave and 
counter-streaming hot electrons with energies of * 1 keV. This interaction occurs in a 
hypothetical region at the geomagnetic equator along the field line of interest in the 
magnetosphere. The interaction region is * 1000 km long and a few degrees of latitude 
wide. Chapter V concludes with a parallel made between whistler amplification and the 
free electron laser (FEL). This discussion includes the idea of a "tapered" interaction and 
the electron pendulum equation that results from a homogeneous geomagnetic field. 
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In Chapter VI, the equations of motion for each hot electron are solved 
simultaneously with the VLF wave equation using computer simulations. These equations 
are numerically integrated given an appropriate set of initial conditions. Chapter VI 
includes simulations for a variety of initial conditions illustrating the effects of gain, 
saturation, the trapped-particle instability and the inhomogeneous magnetic field. The 
simulations are able to produce most of the effects seen in the actual experimental data 
with the exception of triggered emissions. Equations that have traditionally been used to 
predict and distinguish the regions of high and low gain, and strong and weak fields for 
FELs are tested on the VLF wave simulations. Finally, the effects of an initial electron 
energy spread are shown for tapered and non-tapered scenarios. 

Chapter VII deviates from the world of geophysics and enters the world of free 
electron lasers. In this chapter, a short pulse effect, known as limit-cycle behavior, is 
predicted to occur in the Continuous Electron Beam Acceleration Facility (CEBAF) 
infrared (IR) FEL. Limit-cycle behavior is the result of a super-radiant process that occurs 
when the electron pulse length is comparable to the slippage distance. This process causes 
the optical power to oscillate even though all operational parameters are held constant. 
This chapter concludes with a comparison between whistler interactions and free electron 
lasers. 

In Chapter VIII, suggestions for improving the VLF wave simulations are made. The 
major area of improvement includes the incorporation of more realistic distributions 
describing the hot electron’s energy and pitch angle. 
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II. MOTIVATION FOR THE STUDY OF WHISTLERS 
AND FREE ELECTRON LASERS 



A. NAVY VLF COMMUNICATION AND NAVIGATION SYSTEMS 



"Whistlers" are very low frequency (VLF) radio signals in the audio-frequency range 
that resemble a "whistle" [Ref. 2]. Whistlers are characterized by their large sweep in 
frequency over a short time period. A typical whistler can sweep 500 Hz in one second. 
Whistlers can be heard over most audio frequency radio equipment and even long runs of 
telephone line. They are typically initiated by strokes of lightning and propagate in both 
the earth's ionosphere and magnetosphere. The density of whistler signals is dependent 
on the time of day, the earth's magnetic activity, location on earth and the weather. 
Whistlers are more common at mid-latitudes and because of poor propagation conditions 
are not observed at the equator nor the poles. The waveform of a whistler is sketched in 
Figure 2- la. Also shown in Figure 2.1b, is a simple sketch of the usual way whistlers are 
displayed with frequency along the vertical axis and time along the horizontal axis. A 
complete description of all known whistler types is found in Appendix A. [Ref. 2] 



Impulse 



Whistler 
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b 



6000 




4000 



2000 



0 



Figure 2-1. Idealized waveform (a) and spectrum of a whistler (b). [Ref 2] 
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Because the Navy depends heavily on VLF signals, the study of whistlers is 
warranted. The Navy uses VLF radio signals as its primary means for long range 
communication with submarines. Because of the long wavelengths involved, VLF signals 
can travel large distances with low attenuation in the earth-ionospheric waveguide and can 
penetrate sea water to a depth of 10 to 30 feet. Therefore, VLF is used extensively as the 
primarily means of transmitting everything from the news to nuclear release orders to 
submarines. There are many VLF transmitting stations in the United States. They are 
listed in Table 1-1. [Ref. 3] 



TABLE 1-1. VLF TRANSMITTERS. 



Call Sign 


Transmitter 


Frequency 


Latitude 


Longitude 


NSS 


USN Maryland 


21.4 kHz 


39° N 


76° W 


NPM 


USN Hawaii 


23.4 kHz 


21° N 


158° W 


NAU 


USN Puerto 
Rico 


28.5 kHz 


18° N 


67° W 


NAA 


USN Maine 


24.0 kHz 


45° N 


67° W 


NLK 


USN 

Washington 


24.8 kHz 


48° N 


122° W 



Not only are VLF signals used for communications, they are also used extensively for 
long range navigation systems. Both Loran and Omega are navigation systems that 
operate in the VLF, medium frequency (MF), and low frequency (LF) bands. Loran-A, 
originally devised during World War II, operates at 1850-1950 kHz with a range of 
approximately 450 to 800 nautical miles (nm) by day and up to 1400 nm at night [Ref. 4]. 
Loran-C operates at 90 to 110 kHz so that both ground wave and sky wave modes are 
possible [Ref. 5], The ground wave mode has a range of approximately 1200 nm and the 
sky wave mode with two reflections can be received up to 4000 nm. 

Omega operates at 10 to 14 kHz. Since this system operates at very low frequencies, 
it has an effective range of 5000 to 6000 nm and can penetrate sea water so that 
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submarines can use it at periscope depth. Atomic clocks are required to synchronize the 
transmitted and received signals. Only eight stations are necessary to cover the entire 
globe, with three stations being accessible at any location. [Ref. 4] 

B. EFFECTS OF WHISTLERS ON VLF SIGNALS 

As stated earlier whistlers are initiated by lightning discharges known as "spherics" 
which are impulsive in nature. Spherics themselves constitute a major impediment to VLF 
navigation and communication system as a source of intensive static or noise. Suppose 
that a transmitter is sending out a VLF signal with constant amplitude, and a receiver, 
some distance away, is receiving the signal. Typically, the wave energy will propagate in 
the earth-ionosphere waveguide, interacting with the two boundaries in a manner that is 
dependent on the electrical properties of the two boundaries. This is illustrated in Figure 
2-2 [Ref. 6]. Because the electrical properties, such as the conductivity and the relative 
permittivity remain essentially constant over many seconds or minutes, the signal strength 
remains essentially constant and may be synchronized and received. 



tower Ionosphere ( (1 ) 




transmitter receiver 
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Figure 2-2. Schematic representation of a VLF signal (top) propagating under ambient 
conditions with the relative signal strength shown (bottom) [Ref. 6]. 
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When sudden (< Is) or impulsive changes in one or the other waveguide boundaries 
occur, the amplitude and the phase of the signal beyond the disturbance also suddenly 
changes. This can be a significant problem when binary data is being transmitted. A 
sudden change in the lower boundary (i.e. the ground or water) is not as likely, and will 
not be considered here. A sudden change in the upper boundary occurs frequently, and is 
caused by solar flares or by energetic electrons being dislodged from the radiation belts. 
When a whistler or a manmade VLF signal is launched, it can enter the magnetosphere via 
a ducted mode discussed in Chapter III and interact with gyro-resonant electrons. During 
this interaction, the electrons give up energy to the VLF wave and amplify it. If sufficient 
energy is transferred, the electron cannot sustain its original orbit, and is precipitated out 
of the magnetosphere into the ionosphere. Figure 2-3 (top) schematically illustrates this 
type of disturbance. The disturbance grows in size as the precipitated electrons cause 
secondary ionizations where they slow down. Figure 2-3 (bottom) shows a typical VLF 
signature disturbed by a lightning-induced ionospheric disturbance. 





Figure 2-3. Schematic diagram of a perturbed ionosphere (top) and a typical VLF 
signature after a lightning-induced disturbance (bottom) [Ref. 6]. 
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The disturbed region recovers back to its ambient state after approximately 100 s. This 
characteristic signature is commonly referred to as a "Trimpi event", consisting of a quick 
(<1 s) perturbation followed by a relatively slow (*100 s) recovery. 

Therefore, whistlers can affect Naval VLF signals directly and indirectly. Directly, 
they can initiate and sweep over transmitted frequencies resulting in low signal-to-noise 
ratios at the receiving station. Since they are VLF also, they can travel large distances and 
cause communication problems over a large area. Indirectly, whistlers can enter the 
magnetosphere where they can interact with electrons in the Van Allen Belt causing them 
to precipitate. This precipitation can cause Trimpi events which form "pimples" on the 
upper ionosphere waveguide boundary that can disrupt communications for many seconds. 
This thesis concentrates on the magnetospheric interaction between electrons and 
manmade VLF signals, which are used to simulate whistlers. 

C. SHIPBOARD USE OF THE FREE ELECTRON LASER 

With the current threat of third world arms proliferation and the uncertainty 
surrounding the breakup of the Soviet Union, many new countries are achieving the 
capability of launching cruise and ballistic missiles capable of reaching the Unites States 
and its allies. Sophisticated point defense systems such as the Patriot missile (ground- 
based) or the Phalanx Close-In Weapon System (ship based) are not good enough to 
defend against new and improved tactical cruise missiles that are capable of 10 g terminal 
jinking maneuvers. These new missiles could be defended against with a defensive missile 
of roughly 3 times the maneuverability or 30 g's; however, this same technology would 
eventually be applied to the offensive weapon. The escalating "see-saw" between 
offensive and defensive kinetic energy missiles is not practical nor achievable. Therefore, 
the need for a speed of light weapon is obvious. The question becomes, "What type of 
speed of light weapon is best suited to fit the needs of the Navy?" 
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The answer is, of course, a laser. After, the Strategic Defense Initiative (SDI) started 
in the early 1980’s, both neutral and charged particle beams were proven to be inefficient 
methods of transporting energy through the atmosphere. This left the laser as the only 
viable option. Currently, two types of laser systems are under development that are 
capable of generating the necessary power for defensive purposes. They are the mid- 
infrared advanced chemical laser (MIRACL) and the free electron laser (FEL). High 
energy lasers (> 1 MW) can effectively counter high speed and quick reaction targets. A 
missile traveling Mach 2 at a range of 5 km would only close 5.7 mm during the time it 
takes the light to reach it. Typical dwell times required to destroy an incoming missile are 
on the order of 2 seconds. This has obvious advantages over the tracking and homing 
limitations of conventional counter-defense missiles. 

The Mid-Infrared Advanced Chemical Laser (MIRACL) at White Sands Missile 
Range is a deuterium fluoride (DF) chemical reaction laser. The laser output, at 3.8 pm, is 
rated for powers greater than 1 MW. The laser bums NF 2 and D 2 to produce excited 
DF*. The reaction takes place in a resonant optical cavity that couples the coherent 
radiation as the DF* de-excites. The DF gases are exhausted from the cavity using steam 
ejectors. This by-product gas is acidic and highly toxic. Successive testing of the 
MIRACL system has been conducted against supersonic drones. 

On the other hand, the FEL is another high power laser design with many advantages 
over the chemical laser. The FEL uses a relativistic beam of high energy electrons as the 
lasing source [Ref. 7]. Unlike conventional lasers or chemical lasers, a FEL is not a device 
based on atomic transitions. The electrons radiate as they "wiggle" through a periodic 
magnetic field. Therefore, there is no limit to the magazine depth since the supply of 
electrons is endless. It will be shown in Chapter VII, that the wavelength of a FEL is 
proportional to the squared inverse of the electron beam energy. Therefore, the FEL's 
output is tunable over many wavelengths just by adjusting the electron beam energy. 
FELs have demonstrated output wavelengths from the infrared (IR) to the ultraviolet 
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(UV). The advantages of having a tunable laser are numerous. The optimal wavelength 
can be selected so that atmospheric degradation is minimized, and missile skin damage is 
maximized. Finally, compared to a chemical laser, the FEL is a clean weapon. Electrons 
are much more manageable than a hot, acidic DF gas. Therefore, shipboard personnel are 
at less risk while the FEL is operating. 

There are, however, some problems that the FEL must overcome before it could be 
placed onboard a ship. Even though FELs have high peak powers (GW) [Ref. 8], none 
have demonstrated the high average power required for a ship-based weapon. Second, the 
FEL and support equipment require a large amount of space. Therefore, it would be 
impossible to place the system on an existing surface combatant. Finally, to operate a FEL 
with an output of 2 MW and a wallplug efficiency of 20%, the system would require 10 
MW of electrical power. This load would be too large for the onboard turbine generators. 
Therefore, new electrical distribution equipment would have to be designed. However, all 
electrical-drive ship designs, recently being considered, could easily transfer electrical 
propulsion power over to the FEL. Even with these drawbacks, the FEL is a young and 
promising technology that deserves the Navy's attention. 
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III. VLF WAVE-INJECTION EXPERIMENTS FROM 
SIPLE STATION, ANTARCTICA 



A. THE MAGNETOSPHERE 

In order to understand how radio waves and whistlers are amplified, it is imperative to 
accurately describe the environment in which they propagate. The terms and illustrations 
developed in this chapter are important in order to understand the experimental and 
theoretical concepts developed later. Throughout this thesis, the earth's geomagnetic field 
is assumed to equal the field produced by a perfect dipole magnet. This of course is not 
exactly true because of the solar wind. The magnetosphere is broken into two regions 
separated by a boundary known as the magnetopause. The geomagnetic field lines all 
close inside the magnetopause. This boundary separates the region where the earth's 
magnetic field is closely approximated by a dipole from the region where the action of the 
solar wind strongly distorts the earth's field. The region outside the magnetopause is 
known as the magnetosheath. A bow shock wave forms on the sun side of the earth as the 
solar wind causes the magnetosphere to be compressed. As the solar wind passes over 
and around the magnetosphere, a "tail" is formed which has been observed to stretch out 
over hundreds of earth radii. Figure 3-1 illustrates the magnetosphere components and 
how the solar wind distorts the shape. [Ref. 9] 

By concentrating on regions inside the magnetosphere, the dipole magnetic field 
approximation is applicable and accurate. This region changes size depending on the time 
of day, time of year, solar activity, and particle injection events. This region includes the 
plasmasphere which is made up of a "cold" background plasma that determines the bulk 
properties of the media and a "hot" plasma which makes up the earth's Van Allen radiation 
belts. The cold background plasma consists of electrons and protons ranging in energy 
from approximately 0.1 to 1 eV. Their concentration can be described by diffusive 
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equilibrium models with densities ranging from about 10,000 particles/cm 3 at ionosphere 
heights («1000 km) to about 200-400 particles/cm 3 at the plasmasphere boundary. The 
plasmasphere boundary, called the plasmapause, was discovered by using naturally 
occurring whistler waves and is characterized by a significant decrease in cold plasma 
density [Ref. 10]. 




Figure 3-1. A schematic diagram of the earth's magnetic field demonstrating the 
components of the magnetosphere. 

The hot plasma includes electrons and protons having energies ranging from about 1 
keV to 100 MeV which are trapped in helically precessed orbits illustrated in Figure 3-2. 
These energetic particles follow the magnetic field lines and bounce between magnetic 
conjugate points. As the particles move away from the equator their kinetic energy is 
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transferred from parallel velocity components to perpendicular velocity components. A 
convenient parameter to monitor during a particle's trajectory is its pitch angle, 



a = tan ‘ M- 



(3.1) 



where a is the angle between the particle’s perpendicular velocity, V L , and its parallel 
velocity, [Ref. 11]. Therefore, the particle’s pitch angle increases as it moves away 
from the equator until it equals 90° , and then bounces and heads the other way. The 
"bounce latitude" for a particular particle is determined by knowing its pitch angle at any 
particular latitude on its path. The particles also precess around the earth because the 
magnetic field is non-uniform, since it gets stronger closer to the earth. Therefore, the 
particles spiral in non-symmetric helical paths that "walk" or "drift" around the earth. 




-5 x 5-5 x 5 

Figure 3-2. Example of one electron's bounce (right) and drift motion (left) in a dipole 
field where distances are in earth radii. 



A detailed example demonstrating the formulation of these dynamics can be found in 
Appendix B. Typical bounce periods for 1.5 keV electrons range from 4 to 7 seconds and 
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typical drift periods range from 40 to 60 hours. Figure 3-3 plots the bounce period and 
the drift period for a 1.5 keV electron versus its equatorial pitch angle. 




Figure 3-3. Plot of a 1.5 keV electron's bounce period and drift period versus the 
electron's initial equatorial pitch angle at L = 5.1. 

If an electron has a sufficiently low pitch angle at the equator, it will travel far enough 
down its respective field line to enter the ionosphere. Once the electron has entered the 
ionosphere, it can easily interact with many positively charged ions where it can chemically 
bond or lose enough energy from collisions that it will not bounce and return along a field 
line. These electrons are lost from the Van Allen belts due to this depletion mechanism 
known as "precipitation". This limiting equatorial pitch angle subtends a solid angle 
known as the "loss cone" within which there are no electrons. For a given equatorial pitch 
angle, an electron will travel a certain distance into the ionosphere. Figure 3-4 plots the 
equatorial pitch angle an electron must have for a given /--shell to reach an altitude of h = 
1000 km. This "loss height" is somewhat arbitrary but agrees with experimental 
observation [Ref. 11]. The equatorial loss cone pitch angle is insensitive to changes in h 
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as Figure 3-4 also illustrates. The L-shell, or the Mcllwain I-parameter, is a convenient 
way to describe a given magnetic field line in terms of the radial distance at which the field 
line crosses the geomagnetic equator. 



L 



Os. 

K 



(3.2) 



where R eq is the radial distance to the geomagnetic equator and R 0 is the mean radius of 
the earth (6370 km). [Ref. 1 1] The details of Figure 3-4 formulation are also presented in 
Appendix B. 




Figure 3-4. Plot of the electron equatorial loss cone pitch angle versus L- shell where h is 
the ionospheric loss height. 



B. SIPLE STATION EQUIPMENT AND ANTENNA DESCRIPTION 

Very-low-frequency (VLF) transmissions began in Antarctica in 1965 with a 33.6- 
kilometer (km) dipole laid on the ice at Byrd Station (68.2°S, 1=7.25). This facility was 
used jointly by Stanford University and the University of Washington [Ref. 12]. However, 
the geomagnetic latitude of Byrd Station was too high for receiving ducted transmissions 
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in the conjugate hemisphere [Ref. 13]. Accordingly, a lower latitude site (60.8°S, L- 4.2) 
was selected where a 21 -1cm horizontal dipole antenna and a VLF transmitter were placed 
in operation in 1972 at Siple Station. The receiving and recording station was located at 
Lake Missiniti, Canada. This antenna, with a resonant frequency of 5.1 kHz, was used to 
diagnose the magnetosphere for eleven years. In 1983, a new 42-km horizontal dipole 
antenna with a one-half resonant frequency of 2.49 kHz was put into operation. This 
longer, lower wavelength antenna had its resonant frequency in the center of the 1 -4 kHz 
band that had been found to be ideal for wave injection experiments [Ref. 14]. 

In January 1983, current and voltage measurements for the 42-km antenna were made 
for the second resonance frequency of 7.75 kHz. The current distribution throughout the 
antenna allows the ideal transmitted power to be computed assuming the earth is a lossy 
flat plane. The vector potential for the horizontal dipole antenna is computed using, 



where A x is the complex vector potential amplitude in the x-direction, k is the wave 
number, h is the antenna height above ground, R h is the reflection coefficient, and 0 is the 
polar angle in spherical coordinates [Ref. 15]. The complex vector potential amplitude is 
given by [Ref. 15] 



where x is a unit vector in the x-direction, |i is the permeability of the propagation 
medium, r is the far-field radial distance from the center of the antenna, and I(x) is the 
current distribution along the antenna. Once the vector potential is calculated, the electric, 
£, and magnetic, H , fields are obtained and expressed in spherical coordinates as [Ref. 



A = A x (e ,khca<6) + /?,(e)e- , ‘' ,cos(0) ) 



(3.3) 




-ikr L/2 



(3.4) 



-L/2 



15] 



£ = -/© (^0 + ^40 



(3.5) 



and 
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H = -rxE (36) 

n 

where co is the wave frequency and r| is the intrinsic impedance of the medium. The 
Poynting vector follows from the cross product of the real parts of the electric and 
magnetic fields. In this case, the magnitude of the time-averaged Poynting vector, |s|, 
simplifies to 

(3.7) 

where the * denotes the complex conjugate. The Poynting vector magnitude is integrated 
over the upper hemisphere of the antenna, 

P= JJ|S|^ (3.8) 

upper 

hemisphere 

to yield the total power, P, emitted by the antenna. Table 3-1 below gives the 
characteristics of the Siple Station 42-km antenna. [Ref. 12] The radiated powers are 
calculated numerically using Table 3-1 and Equations (3.3) through (3.8) [Appendix C]. 
The radiated powers in Table 3-1 are consistent with the values reported by the STAR 
Laboratory at Stanford University. 



TABLE 3-1. PROPERTIES OF THE 42-km ANTENNA AT SIPLE STATION, 
ANTARCTICA 



Frequency 

(kHz) 


Impedance 

(n) 


Attenuation 

Constant 

(Np/m) 


Phase 

Constant 

(rad/m) 


Phase 

Velocity 

(m/s) 


Radiated 

Power 

(kW) 


2.49 


75 


1.4E-5 


7.5E-4 


2.09E8 


* 1.7 


7.75 


415 


1.8E-5 


2.2E-4 


2.17E8 


*0.31 


13.40 


770 


4.5E-5 


3.7E-4 


2.20E8 


*0.17 
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C. DUCTED VLF WAVE PROPAGATION THEORY 

In order for whistler waves or VLF signals to interact with counter-streaming 
electrons in the magnetosphere, they must follow the same magnetic field lines for an 
extended distance. Typically when a RF wave is launched, either naturally or manmade, 
some of the signal is trapped in the earth-ionosphere waveguide while the rest escapes into 
space depending on the local index of refraction, incidence angle and many other factors. 
However, measurements made at conjugate magnetic receiving stations, show that RF 
waves can travel along "ducts" which form around magnetospheric field lines [Ref. 2]. 
Duct propagation was primarily determined first by the delay time measurements made at 
the receiving station since the launch time was known. 

Electromagnetic waves bend in response to changes in the index of refraction and the 
index of refraction changes as the local electron density fluctuates. Therefore, the 
anisotropy of the ionosphere and magnetosphere at very low frequencies can act to guide 
the path of propagation approximately in the direction of the static magnetic field. 
Evidence has shown that for non-ducted VLF signals, the wave-normal angle becomes too 
large and exceeds the critical angle for total internal reflection. In that case only a weak 
evanescent signal would be received on the ground. [Ref. 2] 

To discuss the theory of VLF wave trapping, a ray tracing duct model is used. The 
duct is formed by enhancements in the index of refraction caused by the electron density 
profile. A useful parameter for studying ducting is the ratio of the whistler wave 
frequency to the electron gyro-frequency, A - col Q.. A special condition occurs when the 
wave frequency is one-half the electron gyro-frequency, A = 0.5. The phase and group 
velocities are equal and the entire character of propagation changes. Most whistler and 
triggered emission phenomena occur near this frequency for reasons not well understood. 
Figure 3-5 illustrates the ducting ray traces for "crest trapping" and "trough trapping." In 



18 



these two figures, N(0) is the electron density on the axis of the irregularity and N(P) is 
the electron density at the outermost excursion of the ray path from the axis. The local 
magnetic field direction is given by H 0 . [Ref. 2] 





Figure 3-5. Ray trace in a field aligned irregularity. Crest trapping, 0 < A < 0.5 (left) and 
trough trapping, 0.5 < A < 1.0 (right). 

The need for ducted propagation was based on three observations of ground-based 
whistlers. (1) A whistler which is excited by a lightning flash appears as a series of 
discrete frequencies with shapes that are independent of the location of the lightning 
source and location of the receiver. This result establishes that the paths of propagation 
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are fixed in the magnetosphere and ionosphere and are a characteristic of the medium 
only. (2) The second observation is closely related to the first and is based on whistler 
echoes. Frequently, a whistler is observed to echo and the separation time of the echoes 
for a given frequency is constant. This suggests that the signal follows the same path each 
time it traverses the magnetosphere. (3) Finally, if a fully developed, well-defined nose 
whistler is carefully analyzed, it fits, within experimental error, the theoretical predictions 
based on purely longitudinal propagation. This close agreement between the shape of 
each trace and the predictions of the purely longitudinal theory requires that the paths be 
field aligned. [Ref. 2] 

D. SIPLE STATION TRANSMITTED FREQUENCY FORMAT 

On January 23-24, 1988 a diagnostic frequency format lasting one minute was 
transmitted from Siple Station every five minutes. This format consisted of a series of 
fixed frequency pulses, frequency ramps, and parabolas between 1.9 kHz and 2.9 kHz. 
The format was centered at the resonance frequency of 2.45 kFiz for the 42-km antenna. 
In addition to the frequency variations, transmitted power was either held constant or 
ramped with a constant slope for some of the signals. Finally, in all cases but two, the 
signals were right-hand circularly polarized and in those two special cases, one signal was 
left-handed and the other was linearly polarized. Each signal was designed with a 
particular experiment in mind in order to probe the different characteristics of the 
magnetosphere. In the next chapter, a detailed study of the macroscopic properties of 
each signal is presented. From this analysis, conclusions are drawn about the properties of 
the magnetospheric duct these signals propagated. Figure 3-6 illustrates the transmitted 
format with frequency (kHz) versus time (sec) during transmission and Table 3-2 gives the 
details about each signal. 



20 




Figure 3-6. Example of the HR241 Siple Station diagnostic format. The center frequency 

is 2.45 kHz. 



TABLE 3-2. HR241 SIPLE TRANSMITTER FORMAT AND SIGNAL 

CHARACTERISTICS. 



Signal 


Frequency 

(kHz) 


Power Level 


Polarization 


Signal 

Length (sec) 


Time 

Transmitted 


la 


2.4±.5 ramp 


10 dB/sec 


Right Hand 


1 


0:00 


lb 


2.4 constant 


for 1/2 sec 




1 


0.01 


2a 


2.4±.5 ramp 


10 dB/sec 


Left Hand 


1 


0:02 


2b 


2.4 constant 


for 1/2 sec 




1 


0:03 


3a 


2.4±.5 ramp 


10 dB/sec 


Linear 


1 


0:04 


3b 


2.4 constant 


for 1/2 sec 




1 


0:05 


4 


2.4 ±.25 
ramp 


constant 


Right Hand 


4 


0:13 


5 


2.4 constant 


constant 


Right Hand 


2 


0:17 


6 


2.4 constant 


constant 


Right Hand 


200 msec 


0:38.5 


7 


2.4 constant 


10 dB/sec 
for 1 sec 


Right Hand 


2 


0:41 
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IV. ANALYSIS OF SIPLE STATION DATA FROM 23-24 
JANUARY 1988 



A. WHY WAS THIS DATA SET CHOSEN? 

Between 1700 UT (universal time) on 23 January to 0300 UT on 24 January 1988, 
the HR241 format described in Chapter III was transmitted from Siple Station, Antarctica. 
This data set was selected from thousands of hours of transmitted data because of many 
desirable characteristics. First of all, most signals that are transmitted from such a long 
non-directional antenna excite many ducts. This occurs because poor directivity allows 
many transmission angles to receive significant RF power. Also, the dispersive ionosphere 
can bend and scatter signals into many magnetospheric ducts. Therefore, the receiving 
Station receives signals that may have propagated through many ducts with different delay 
times. This is called ’’multipath” propagation as illustrated in Figure 4-1. In this figure, 
the top panel shows a whistler that propagated over three distinct paths with three distinct 
reception time delays. This is usually not desired because it is difficult to analyze what 
happen in just one duct. The bottom panel shows a pure, single-path whistler with its 
echo approximately one second later. This data set showed little multipath propagation 
over the ten hour period. 



*.p, *. a,. B»., B„ 9«. 
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Figure 4- 1 . Example of a multipath whistler (top) versus a single-path whistler (bottom). 
Frequency is along the vertical axis and time is along the horizontal axis [Ref. 2]. 
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To compare data that spans over a ten hour period, it is important that all signals 
traverse the same duct and stay in the same duct. This is determined by first observing 
single-path propagation and second, measuring the signal delay time. The signal delay 
time is the time difference between transmission and reception. Consistently throughout 
the ten hour period, the delay time of each signal was four seconds. This was true when 
different signals were compared to each other within the one minute transmission period 
and when the same signals were compared to each other after each five minute off period. 

These data were also interesting because of some atypical time variations in the type 
and number of signals received. Moreover, there was an extremely low density of 
triggered emission activity observed during this ten hour period. It is important and 
helpful to graphically study a few examples: (1) There were times during a one minute 
transmission that only frequency ramped signals were received and a few cases when only 
constant frequency signals were received as illustrated in Figure 4-2 and 4-3. 




Times Past 17:05 (UT) on 23-24 Jan. 1988 



Figure 4-2. Times during the one minute transmission period when only ramped 
frequency signals were received. 
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Figure 4-3. Times during the one minute transmission period when only constant 
frequency signals were received. 



At the present, there is no theory or adequate explanation to describe why one type of 
frequency profile was preferentially amplified and received. When this rare behavior is 
observed, it is more likely that a frequency ramped signal will be seen. This occurs 
because the signal is sweeping over a large frequency band which increases the probability 
that the VLF wave will find a suitable resonance condition with the interacting hot 
electrons. (2) Within a one minute transmission period, there was a total of 13 signals 
transmitted. The total number of signals received at Lake Missiniti, Canada out of 13 
transmitted varied significantly over the ten hour period. The percentage of signals 
received in a one minute period over the ten hour time scale is shown in Figure 4-4. 



24 




Figure 4-4. Percentage of the total signals received in a one minute transmission period. 



Using Figure 4-4 as a guide, an hour and a half window starting at 00:15 (UT) was chosen 
to be studied in more detail because of the large percentage of transmitted signals received 
during this period. Even over a one minute time scale not all signals of the same type 
were received. Therefore, the mechanism for amplification can change abruptly during the 
one minute transmission period. However, Figures 4-2 and 4-3 depict that this was 
usually not the case. One has to factor in that these signals may or may not be at the same 
initial power (see Table 3-2) when analyzing Figure 4-4. (3) Finally, it is typical for Siple 
data to be plagued with many triggered emission signals sprouting off transmitted signals. 
"A triggered emission is any emission that appears to have been initiated, or triggered, by 
another event or interaction such as whistlers, discrete emissions, or signals from VLF 
transmitters [Appendix A]." Emissions can easily sweep fifty percent of their initiation 
frequency and usually arrive in multiple ducts. Figure 4-5 demonstrates a transmitted 
signal with rising and falling triggered emissions emanating from it. However, the data 
received during the discussed ten hour period had a very low density of triggered 
emissions and all but a very few cases the emissions were single-path. Again, this made 
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analyzing this data set easier. The triggered emissions usually occurred at the termination 
of the transmitted signal and demonstrated very pure frequency sweeps with no sidebands. 
These rare characteristics demanded a closer and in-depth study of the data set from 23-24 
January, 1988. 




Figure 4-5. An example of rising and falling triggered emissions [Ref. 2]. 

B. RECEIVED FREQUENCY AND AMPLITUDE INFORMATION. 

The receiving station recorded all of the data on magnetic analog tapes that are stored 
at Stanford University. This raw data was processed with a 100-Hz band pass filter 
around the center frequency of 2.45 kHz and then plotted on amplitude versus time 
spectrograms as demonstrated with Figure 4-6 (bottom). The amplitude data was 
obtained by taking the fast-Fourier transform of the data and applying a low pass filter in 
order to eliminate high frequency '’atmospherics" that are always present during data 
collection. After removing the high frequency noise components from the data, the 
inverse Fourier transform was performed. These atmospherics, or "spherics” for short, are 
strong impulses produced by stokes of lightning and are what typically initiate whistlers. 
Examples of spherics are seen as dark horizontal striping in Figure 4-5. In Figure 4-6, the 
frequency and amplitude information are shown along the same time-scale axis where the 
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large tick marks are one second apart The center of the dark frequency band is at 2.45 
kHz and the maximum width is approximately 200 Hz and the triggered emission falls 
approximately 600 Hz. The amplitude of the signal saturates at approximately -10 dB 
where -8 dB corresponds to 300 jiV/m. Saturation is a direct result of non-linear wave- 
particle interaction. 




Figure 4-6. Example of Signal #5 (constant frequency and power) illustrating growth, 
saturation, sidebands and a triggered emission. Frequency vs. time (top) and Amplitude 
vs. time (bottom) are against the same time scale where the tick marks are one second 

apart. 
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1. Exponential Growth 

As the VLF wave travels along a magnetic field line it may be amplified by counter 
streaming energetic electrons that spiral around the same field line. The details of this 
interaction are presented in the next chapter. The exponential growth rate for this 
example is 35 dB/sec and is typical of the signals seen during the ten hour period. The 
exponential growth continues until the signal saturates. Notice from Figure 4-6 that the 
frequency (top panel) remains relatively narrow during the initial growth period, and then 
broadens as sidebands are formed during saturation. The growth rates varied significantly 
over the ten hour period as demonstrated by Figure 4-7 which displays a 1.5 hour 
window. This figure demonstrates large changes in growth rate over a five minute time 
scale for Signal #5. 
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Figure 4-7. Plot of growth rate (dB/sec) versus time for Signal #5. 



To investigate changes over a shorter time period, it is necessary to look at Signal #7 
which was transmitted 24 seconds later. Figure 4-8 plots the ratio of Signal #7 growth 
rate to Signal #5. Signal #5 and Signal #7 were chosen as pairs because they are the most 
similar signals transmitted. This figure illustrates that the growth rates for these two 
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signals are the same within the uncertainty of the measurement given by the error bars. 
The average growth rate ratio was 0.95 ±0.19. It should also be noted that the growth 
rates were difficult to measure precisely and contain an average uncertainty of ± 5 dB/sec. 
From Table 3-2, Signal #5 was transmitted with a constant power level; however, Signal 
#7 was ramped in power at 10 dB/sec. Therefore, one may ask how can these two signals 
be compared if one has a transmitted power ramp in its profile? The answer to this will be 
discussed fully later in this chapter. 




Figure 4-8. Ratio of Signal #7 (corrected) growth rate to Signal #5 growth rate during a 

1.5 hour window. 



2. Saturated Power Level 

Typical to any system with negative feedback, an amplified signal reaches a steady- 
state power level. Therefore, this section will discuss some of the properties of the 
observed saturated levels and the next chapter will develop the wave-particle interaction 
and feedback mechanism in detail. Once again the saturated power level for Signal #5 
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varied significantly over a five minute time scale as illustrated by Figure 4-9. All of these 
saturation levels are normalized to what is thought to be the lowest unamplified power 
level during this 1.5 hour window. The calculation of the unamplified power level is not a 
trivial matter and will be explained fully later in this chapter. Note the normalized power 
levels, Pf/P 0 in Figure 4-9 can also be thought of as gain levels, 



G = 



Pr~P. 






(4.1) 



where G is the signal gain, P^-and P 0 are the final power and initial unamplified power level 
respectively. This is an excellent approximation when Py » P 0 . 
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Figure 4-9. Chart of normalized power levels for Signal #5. 



It is easy to see that gains as large as 35 dB are possible and that gains on the order of 30 
dB occur frequently. Large gains like this are one of the main reasons these signals are 
studied. If the mechanism of this gain process can be fully understood, similar principles 
could be applied to free electron laser theory and design. 
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It is interesting to see how the saturated power levels change on the 30 second time 
scale. Just as was done for growth rate above, the ratio of the saturated power level for 
Signal #7 to that for Signal #5 is plotted in Figure 4-10. 




Figure 4-10. Chart of the ratio of saturated power level for Signal #7 to Signal #5. 

When both signals reached a saturated level within the one minute transmission period, 
their signals levels are approximately equal, and thus the ratio is one. Using Chauvenet’s 
criterion [Ref. 16] to reject outlying data in Figure 4-10, the average ratio was 1.02 ± 
0.04. When one of the signals did not grow to saturation for one reason or another, no 
correlation between Signal #5's and Signals #7's saturated signal level was observed, and 
the ratio was not constant nor predictable. In general both signals did saturate and on the 
30 second time scale, their saturated signal levels are nearly equal. Estimating the time 
scale over which key experimental parameters change is useful in understanding the 
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mechanisms that cause the change. For this reason, the time scale question will continue 
to be emphasized throughout this chapter. 

3. Sidebands 

Sidebands produced during the ten hour period had the least variation of all of the 
experimental phenomena observed. Sidebands were only observed if the transmitted 
signal saturated, and were weakly observed in the fine structure of the subsequent 
triggered emission. During the saturated portion of the received signal, up to 8 harmonics 
of the sideband modulation were observed. The sideband frequencies were calculated by 
two methods. First, after the signal saturated, the up and down amplitude peaks, seen in 
the bottom panel of Figure 4-6, were counted and converted into a frequency. Second, 
the Af was read directly off the frequency-time histograms. Both of these methods 
produced the same result within ± 5 Hz. However, the second method is far easier and 
more accurate because the amplitude variations are not simply sinusoidal because they 
contain all of the sideband information. Therefore, it is not obvious which peaks 
correspond to which sideband. The top panel of Figure 4-6 clearly illustrates the sideband 
broadening with Af » 25 Hz. The reader may compare the results with the bottom panel 
of Figure 4-6. To illustrate the invariance of the sideband formation. Figure 4-1 1 plots the 
first sideband of Signal #5 versus time during a 1.5 hour window. The average 
modulation frequency for the first sideband of Signal #5 was 23.2 ± 1.8 Hz. 
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Figure 4-11. Plot of the first sideband frequency Af versus time for Signal #5. 



The formation of sidebands during an amplification process can be a complex subject. 
Sidebands are produced when the main carrier frequency is modulated by a secondary 
interaction process causing the signal to "beat." The modulation or beat frequency (Af) 
can be many orders of magnitude less than the carrier frequency as is the case with the 
Siple experimental data as well as common FM radio signals. A clear cause for the 
sideband production in the Siple data is still under investigation. Possible explanations 
include (1) modulation by the transmitting antenna’s peripheral equipment, (2) modulation 
by Canada's power line grid during reception, and (3) modulation caused by synchrotron 
oscillations due to strong field coupling of the counter-streaming electrons. The last 
possibility will be described in the next chapter and simulations will demonstrate good 
correlation with observed sideband frequencies. 

4. Triggered Emissions 

There are many different types of triggered emissions as already noted in Appendix A. 
During the ten hour observation period, over 90% of the triggered emissions were 
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"termination emissions." Termination emissions are triggered emissions that occur at the 
end of the transmitted pulse. Figure 4-6 is a good representative example of the type of 
termination emissions observed. When the signal is terminated, an emission is launched of 
equal signal strength that first rises in frequency then sharply drops in frequency. The rise 
in frequency is about 5% of the center frequency (2450 Hz); however, the drop in 
frequency is much larger and ranges between 20 to 50% of the center frequency. Once 
again no theory nor simulation has been able to describe triggered emission sweeps of 500 
to 1200 Hz. Similarly, no theoretical mechanism has been formulated that explains why 
triggered emissions occur, and why they occur preferentially on the termination of a 
transmitted signal. One simple macroscopic explanation for the falling emission may be in 
terms of "phase-locked" hot electrons. Once the injected VLF wave train passes by the 
counter-streaming electrons, they are "forced" together in phase due to the influence of 
the wave. The electrons stay in phase as they approach the geomagnetic equator where 
the magnetic field, and therefore the gyro-frequency decreases. The electrons continue to 
radiate in proportion to the local gyro-frequency until they have lost sufficient energy to 
lose their phase-lock, or until they pass through the geomagnetic equator where the 
magnetic field strength gradient changes sign and disrupts the process. The process of 
electrons radiating coherently in the absence of an external electromagnetic field is called 
"super-radiance." There is a free electron laser (FEL) analog to this process which will be 
discussed in Chapter VII; however, the time scales are much different. For the Siple 
Station data, super-radiant coherency would have to be maintained for 0.2 to 0.5 seconds 
to explain the falling triggered emissions in all cases. But for FELs, it has been shown that 
electrons will only maintain their coherency in the absence of light for times on the order 
of nanoseconds. 

It is easy to determine whether or not changes in the electrons' gyro-frequency are 
the same order of magnitude as the falling triggered emission sweeps. Consider electrons 
on the L = 5.1 magnetic field line as is the case for the Siple Station data, and suppose 
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